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Abstract 

We establish a series representation of the Hill discriminant based on the spectral 
parameter power series (SPPS) recently introduced by V. Kravchenko. We also 
show the invariance of the Hill discriminant under a Darboux transformation 
and employing the Mathieu case the feasibility of this type of series for numeri- 
cal calculations of the eigenspectrum. 

Keywords: Hill's discriminant, spectral parameter power series, supersymmet- 
ric partner equation, Mathieu equation. 



1 Introduction 

There is recent strong interest in differential equations with periodic coefficients 
due to their numerous applications in modern material sciences and engineer- 
ing. The main goal in their theoretical framework is the description of the 
spectrum and getting the periodic and quasiperiodic solutions. In the case of 
linear second-order ordinary differential equation, a function of the spectral pa- 
rameter known as Hill's discriminant is the basic quantity containing important 
information both about the spectrum of the differential operator and also about 
the construction of the (quasi)periodic solutions. 

This paper focuses on a representation of the Hill discriminant which is given 
in the form of a power series in the spectral parameter. This representation is 
obtained by using recent results of Kravchenko [H HI 13] • It is worth mentioning 
that since the Hill discriminant for each value of the spectral parameter can 
be obtained from a couple of corresponding linearly independent solutions, any 
available representation for these solutions leads to a representation of Hill's 
discriminant. Nevertheless, in general no easy and practically treatable repre- 
sentation for Hill's discriminant is known unless for the case of some well-studied 
equations such as the Mathieu and the Lame equations. 
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Long ago, Jagerman [3] introduced and studied in detail the so-called car- 
dinal series representation of Hill's discriminant. In the famous book of Mag- 
nus and Winkler [5] the Hill discriminant for Schroedinger type equations is 
expressed as an infinite determinant involving the Fourier coefiicients of the po- 
tential as well as the spectral parameter. The phase-integral method is used 
by Froman [B] to obtain a representation involving a matrix whose entries are 
complicated phase integrals, while Boumenir [7j wrote it in terms of integrals 
derived from the inverse spectral theory. 

In all the aforementioned series representations the spectral parameter enters 
in a quite sophisticated way with all the terms being functions of the spectral 
parameter. Instead, our result herein gives the Hill discriminant in the form of 
a power series in the spectral parameter with the series coefficients independent 
of it and calculated only once, i.e., for a single value of the spectral parameter. 
Moreover, its practical implementation is easy and as an illustration we show in 
the present work the numerical results obtained for the Mathieu equation. 

We also show that the SPPS representation gives additional information not 
only concerning the original equation but also on its Darboux-related partners. 
As a corollary, we prove the invariance of the Hill discriminant under the Dar- 
boux transformation under some additional conditions. 

2 Hill's type equations 

The Sturm-Liouville differential equation 

L [f{x, A)] = -{p{x)nx, A))' + q{x)f{x, A) = A/(x, A) (1) 

with T-periodic coefficients p{x) and q{x) and real parameter A is known as 
of Hill type. We first recall some necessary definitions and basic properties 
associated with the Eq. ((!]) from the Floquet (Bloch) theory. For more details 
see, e.g., [5j[8]. In what follows we assume that p{x) > 0, p'{x) and q{x) are 
continuous bounded functions. 

For each A there exists a fundamental system of solutions, i.e., two linearly 
independent solutions of dJ) /i(a;. A) and /2(x, A) which satisfy the initial con- 
ditions 

/i(0,A) = l, /((0,A)=0, /2(0,A)=0, /^(0,A) = 1. (2) 

Then the Hill discriminant associated with Eq. ^ is defined as a function of A 
as follows 

7^(A) =/i(T,A)-l-/^(T,A). 

Employing D(X) one can easily describe the spectrum of the corresponding 
equation. Namely, the values of A for which |I?(A)| < 2 form the allowed bands 
or stability intervals meanwhile the values of A such that |£'(A)| > 2 belong 
to forbidden bands or instability intervals [5]. The band edges (values of A 
such that |-D(A)| = 2) represent the discrete spectrum of the operator, i.e., they 
are the eigenvalues of the operator with periodic {D{\) = 2) or antiperiodic 
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{D{X) = —2) boundary conditions. The eigenvalues A„, n — 0,1,2,... form an 
infinite sequence Aq < Ai ^ A2 < A3..., and an important property of the mini- 
mal eigenvalue Aq is the existence of a corresponding periodic nodeless solution 
fo{x, Ao) [S]. In general solutions of (P) are not of course periodic, and one of the 
important tasks related to Sturm-Liouville equations with periodic coefficients 
is the construction of quasiperiodic solutions. In this paper, we use the match- 
ing procedure from [9] for which the main ingredient is the pair of solutions 
fi{x,X) and f2{x,X) of H]) satisfying conditions Namely, using A) and 
f2{x,X) one obtains the quasiperiodic solutions f±{x + T) = (3±f±{x) as follows 

f±{x, A) = PlF^{x - nT, A), | "J|o^^i|±J^ ■ (3) 

where F±{x, A) are the so-called self-matching solutions, which are the following 
linear combinations F±(x, A) = fi{x, X) + a±f2{x, A) with a± being roots of the 
algebraic equation /2(T, X)a'^ + (/i(r. A) - /^(T, A))a- /{ (T, A) = 0. The Bloch 
factors (3± are a measure of the rate of increase (or decrease) in magnitude of 
the self-matching solutions F±{x, A) when one goes from the left end of the cell 
to the right end, i.e., /3±(A) ~ ^^^(0 a) ' "^^^ values of /3± are directly related 
to the Hill discriminant, P±{X) = ^{D{X) ^ y^D'^{X) — 4), and obviously at the 
band edges f3+ = ~ ±1 for D{X) — ±2, correspondingly. 

3 The SPPS series representation of Hill's dis- 
criminant 

In this section, we will give an efficient representation for the Hill discrimi- 
nant using the method of spectral parameter power series (SPPS) [2 [3], which 
was used in |10[ 111] for studying quantum mechanical models related to one- 
dimensional Dirac systems, also for solving electromagnetic scattering problems 
[12] and for solving inverse problems [13]. The SPPS method offers a proce- 
dure for constructing the solutions fi{x, A) and f2{x, A) of ^ which satisfy the 
initial conditions ^ . This construction is based on the knowledge of one non- 
vanishing particular solution /o(a;,Ao) of ([T]) being bounded on [0,r] together 
with • '^^^'^ '^^ Hill's equation the first eigenvalue Aq of the Eq. (|T]) 

generates nodeless periodic eigenfunction fo{x, Aq). In what follows we initially 
suppose that the value of Aq is known. Note that, it can be obtained by different 
methods including the same SPPS method [3] as we explain in subsection 13.41 

Given Ao , we proceed in three steps in order to obtain the representation of 
Hill's discriminant: 

i) the first one is the construction of a particular nodeless solution 
foix, Ao) of (P which is periodic, i.e., fo{x + T, Ao) = fo{x, Ao), 

ii) the second one is the construction of the fundamental system of 
solutions fi{x,X) and /2(x, A) of ([T]) for all values of the parameter 
A, 
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iii) the final step is getting the representation of Hill's discriminant. 
We detail each of the steps in the following subsections. 



3.1 The nodeless periodic solution 

In order to obtain the nodeless periodic solution fo{x, Xq) of ^ for A = Aq, i.e., 

- {p{x){fo{x)yy + q{x)fo{x) = Ao/o(x) (4) 

we first construct two linearly independent SPPS solutions of (|4]), /o,i(x, Aq) 
and /o. 2(3^5 Aq), which satisfy the initial conditions /o.i(0,Ao) = /o2(0iAo) = 
1 and /o^i(0, Ao) = /o,2(0, Aq) — p. . These solutions are not necessarily 
periodic. The procedure of James jSJ allows one to obtain from /o,i(a;,Ao) 
and fa.2{x, Aq) the Floquet type solutions which degenerate to a single periodic 
solution fo{x, Aq) since Aq represents a band edge. 

The functions fo,i{x, Aq) and /o,2(a;, Aq) can be calculated as follows [1] 



piO) J2 ^o"\ (5) 

odd n—1 



for an odd n 
for an even n 



fo,i{x,Xo)= ^ ^o""* and /o,2(a;,Ao) 



even n—0 

where 



Xt\x) 



xir\x) 



lo ^0" ^^(O^t^^ for an odd n 

lo ^o""'^(0(9(0 - Ao)d^ for an even r 



Now following [3] we obtain the periodic nodeless solution of (|1]) 

/o(a;, Ao) = fo,i{x - nT,\o) + a/o,2(a; " nT, Aq), (6) 

nT <x < {n + l)T 
n = 0,1,2,... . 

where a = ^'^2f°l (tAq)"^ '^"'' since Aq is the band edge eigenvalue [9]. 
3.2 Fundamental system of solutions 

Once having the function fo{x, Aq) the solutions fi{x, A) and f2{x, A) of ([T|) and 
([2]) for all values of the parameter A can be given using the SPPS method once 
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again [T] 

/i(a;,A) = ^So(x,A,Ao)+p(0)/^(0)/o(a;)I]i(x,A,Ao), 

(7) 

/2(x,A) = -p(0)/o(0)/o(a;)Si(x,A,Ao). 
The summations J^q and Zli are the spectral parameter power series 

oo oo 

Eo(x, A, Ao) = ^^'"H2;)(AA)", J:^{x, A, Ao) = ^ x(2"-i)(a;)(AA)"-i , 

n— n— 1 

where AA = A - Ao and the coefficients X<-''^{x), are given by the 

following recursive relations 

= 1, EE 1, 

f r X("-'H0fS{0d^ for an odd n 
X(-\x) = i (8) 

foranevenn 

rjj=^(„-i)(^)_<|_ foranoddn 

^ <^ (9) 

[/Q"X("-i)(e)/o2(0'^C foranevenn . 

One can check by a straightforward calculation that the solutions /i and /2 
fulfill the initial conditions for this the following relations are useful 

(So(x,A,Ao)) = T^X ^^ whereEi(x,A,Ao) = ^x(2"-i)(x)(AAr 

V /x P[X)Jq[X) 

(10) 

and 

(Si(x,A,Ao)): = - "r":^^, where So(x,A,Ao) = ^ x(2")(a;)(AAr . 

(11) 

Having obtained the fundamental system of solutions for any value of A, one can 
apply the construction (jS]) in order to obtain the Bloch solutions which become 
eigenfunctions for A being eigenvalues. 
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3.3 Hill's discriminant in SPPS form 



Now the Hill discriminant D{X) — fi{T, A) + f2{T, A) can be written in a simple 
explicit form. For this we write fi {T, A) and /2 (T, A) in a form of a spectral 
parameter power series using ([7]) and taking into account 

+ (/^(0)/o(T) - /o(0)/^(r))p(0)Ei(r,A, Ao) . 

Finally, taking into account that foix) is a T-periodic function /o(0) = fo{T) 
and writing the explicit expressions for Eo(T, A, Aq) and T,q{T, A, Aq) we obtain 
a representation for Hill's discriminant associated with ([l} 

oo 

D{X) = (l(2n)(y) ^ x(2«)(r)) (AA)". (13) 

n=0 

Thus, only one particular nodeless and periodic solution /o(a;, Aq) of ([J) is needed 
for constructing the associated Hill discriminant. We formulate the result (|13p 
as the following theorem: 

Theorem 1 Let Aq he the lowest eigenvalue of the Sturm- Liouville problem 
for the operator L on the segment [0,T] with periodic boundary conditions and 
/o(x, Aq) be the corresponding eigenfunction. Then the Hill discriminant for {Ip 
has the form US\) where X*-^"^ and X'^"-* are calculated according to ^ and 
(0), and the series converges uniformly on any compact set of values of X. 



To illustrate the formula (|13p we consider a simple example. Let in equation 
(dl q{x) = 0, p{x) = 1. It is easy to see that the associated discriminant is 
D{X) = 2 cos a/at, from where we obtain Aq = and a corresponding non-trivial 
periodic solution is /o(a;) = 1- Now making use of this solution we construct 
the discriminant by means of the formula i^^. The coefficients X(2")(r) and 
X(2n)(T) given by ® and © take the form 

X(2")(r) ^ X(2")(T) = i-irjTr^, = 0, 1, 2, ... . 

(2n). 

The substitution in ^ gives D{X) = 2cos V^T. 



3.4 Construction of the first eigenvalue Aq by the SPPS 
method 

Notice that in the expression ([T^ for D{X) and in all reasonings previous to 
it we do not use the periodicity of the solution /o(x, Aq), therefore and 
the whole procedure for obtaining it are valid for any A* such that there exists 
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a corresponding solution f^,{x,X^) which is bounded on [0,T] together with 
^/(pf*)- Such a solution /, (x, A,) can be obtained in the following way 

/*(a;, A,) = /*4(a;, A*) + if*^2{x, K) (14) 

where /*,i(a;, A,) and /*.2(a;, A,) are given by ([5]) with A, instead of Aq. For more 
details see The pair of the independent solutions fi{x, A) and f2{x, A) of (dJ 
given by (O of course are independent of the choice of the solution /o(a;, Aq), 
hence instead of /o(a;, Aq) in (IT2]) one can take /*(a:, A*) given by (IT4l) . Thus, 
in terms of /,(a;, A*) where A» is essentially arbitrary, I?(A) can be represented 
as a series in powers of (A — A») 

+ if'MMT) - /.(o)/:(r))p(o)x(2"+i)(r)^ (A - a,)". 

Now the band edge Aq required for the formula (|T3)) can be calculated as 
a first zero of the expression D{X) — 2 where D{X) is given by ([T5|) . For the 
numerical purpose it can be useful to know the interval containing Ao- Since q is 
a bounded periodic function, there is a number A which satisfies the inequality 
q{x) > A Vx e R. It is known [S] that D{\) > 2 for aU A < A, therefore the 
lower estimate for Aq is the following 

Ao ^ min q{x). 

The upper bound can be obtained considering the Rayleigh quotient for periodic 
problem [14] 

{u'{x)f + q{x) {u{x)f)dx 

Aq ^ -f ^ , 

Jo ("(2;)) dx 

where u{x) € [0,7"] is periodic with period T. The equality occurs if and 
only if u{x) is an eigenfunction corresponding to Aq. 

4 Hill's discriminant of the supersymmetric (SUSY)- 
r elated equation 

In this section, we consider the SUSY-related equation of Eq. ([T|) and obtain 
its SPPS solutions. These solutions allow us to prove the equality between the 
Hill discriminants of equation ([T]) and its SUSY-related Eq. (IT51) . For various 
aspects of SUSY periodic problems, see [T5l [TBI fT7) . 

The left-hand side of the equation ^ can be factorized in the following way 

m 

L [fix, A)] = (a;) + + $(2;)) f{x), (16) 
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where means the x-derivative, the superpotential ^{x) is defined as fohows 
${a;) = —pi (x) ■ Using this factorization the coefiicient q{x) can be 

expressed as 

q{x)^<i>^{x)~ (^p^x)^x)^ +Xo. 
Introducing the fohowing Darboux transformation 

(pi {x)d, + ^x)) fix, A) = /(x, A), (17) 

one obtains the equation supersymmetrically related to equation ([!]) 
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'fix, A)] = [pi ix)d, + $(x)) (-d^pi ix) + $(.t)) fix, A) = A/>, A), 
which can be written as foUows 

- d.Apix)dJix, A)) + qix)fix, A) - Xfix, A), (18) 
where g(x) is the SUSY partner of the potential qix) given by 

qix) = qix) +2piix)<^'ix) -piix)ipiix))". (19) 

It is worth noting that as "I>(a;) is a T-periodic function, the Darboux transfor- 
mation assures the T-periodicity of qix). 

The pair of linearly independent solutions /i(x, A) and f2ix. A) of can be 
obtained directly from the solutions ([7]) by means of the Darboux transformation 
P7)) . We additionally take the linear combinations in order that the solutions 
/i(a;, A) and /2(x, A) satisfy the initial conditions /i(0, A) = /2(0, A) = 1 and 
/((0,A) = /2(0,A) = 

f I ,^ P^(0)/o(0) ^ , ^ [piix)]'\,^^--fiO) ~ 

jiix,\) — —, Lo(.i;,A)H J LiiXjA), (20) 

pHx)foix) iAX) foiO)pHx)foix) 

f2ix,X) = Si(x,A). (21) 

iAX)foiO)pHx)foix) 

These two solutions allow us to write the expression for Hill's discriminant 
associated to the equation (fTS)) . that is i3(A) = /i(T, A) + f2iT, A). We consider 
first the derivative of f2ix, A) and evaluate it for x = T 

~f'(r All - ^h(n^ \ph^)'\'\^-Tf,iT)+pHT)fl,iT) ~ pH0)foiT) ~ 

f,ix, A)u., - -p^ (0) ^X)7»(^W^^ ' 

Notice that as the functions foix, Xq) and p(x) are T-periodic, i.e., /o(0, Aq) = 
/o(T', Ao) and p(0) — p(T), then obviously, the functions /q(x, Ao),p2 (x) and 
[^2(2;)]' possess the same properties. Therefore we have 



i?(A) = i]o(r, A) + So(r, A) + 



[pHx)]'u^o-m 



(AA)/o(0,Ao)p5(T)/o(r,Ao) 
[pHx)]'U=TfoiT, Xo)+pHT)f(,iT,Xo) 



(AA)/o(0,Ao)p^(r)/2(r, Ao) 



Si(T,A). 
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The substitution $(0) = -p^ (0) -^"["'^"l clearly shows that the expression m 
brackets vanishes. Therefore, we obtain 

oo 

5(A) = I]o(T, A) + So(T, A) = ^ (^X^^"\T) + x(2")(r)) (AA)". 

Comparing with (jl3p we have the identity 

D{X) = D{X). (22) 

Thus, we have proven the following statement; 

Theorem 2 Let Xq be the first eigenvalue of [l^and fo^x, Xq) the corresponding 
T -periodic nodeless eigenfunction. Then the Darboux transformation |J7| ) with 
^(x) = — (x) leads to a SUSY-related Eq. il8\) with the preservation 

of the Hill discriminant, i.e., Eq. h22fl holds. 

From the identity of discriminants ((22|) it is clear that Aq gives rise to a 
nodeless periodic solution fo{x, Aq) of Eq. ([T^. Taking A = Ao in (pn| and (PT|) 
we get this eigenfunction in the form fo{x, Aq) = —i . 

p7(a;)/o(a:,Ao) 

Notice that, the factorization method can be applied to Eq. (|18p with the 
superpotential $i(x) = — (x) {"^^'^"^ . In this case, we obtain the representa- 
tion ^ 

q^^lix) - (p5(a;)$i(x)) +Ao, 

which reduces to the equality ([T9|) if one notices the relationship ^i{x) = 
(p2 (x))' — $(a;). It can be also shown that q = q, where q = q{x)-\-2p2 (x)^[{x)~ 
p^ {x){p^ (x))" is the superpartner potential of q{x). Thus, the Darboux trans- 
formation (fTT]) with the superpotential ^i{x) applied to Eq. (fT8|) does not pro- 
duce a different potential. 

5 Numerical calculation of eigenvalues based on 
the SPPS form of Hill's discriminant 

As is well known, see e.g., [5], the zeros of the functions D{X) ± 2 represent 
eigenvalues of the corresponding operator. In this section, we show that besides 
other possible applications the representation (|13l) gives us an efficient tool for 
the calculation of the discrete spectrum of a periodic Sturm-Liouville operator. 

The first step of the numerical realization of the method consists in calcula- 
tion of the minimal eigenvalue Aq by means of the procedure given in subsection 
13.41 and subsequently in construction of the corresponding nodeless periodic so- 
lution fo{x, Xq) using formula ([6]). The next step of the algorithm is to compute 
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the functions and X^") given by ^ and respectively. This construc- 
tion is based on the eigenfunction /o(a;, Aq). Finally, by truncating the infinite 
series for D{X) we obtain a polynomial in A A 

N 

Dn (A) - J2 (^^'"^ + ^^'"^ (^)) (2^) 

n=0 

TV 

= 2 + ^ (x(2")(T) + X(2")(T)) (AA)". 

n=l 

The roots of the polynomials Dn{\) ± 2 give us eigenvalues corresponding to 
Ec^. ([T]) with periodic and antiperiodic boundary conditions. 

As an example, we consider the Mathieu equation with the following coeffi- 
cients 

p{x) ~ 1, q{x) — 2rcos2a;. 

The algorithm was implemented in Matlab 2006. The recursive integration 
required for the construction of ^q"-*, Xq"', X^"^ and X^"^ was done by rep- 
resenting the integrand through a cubic spline using the s'pa'pi routine with a 
division of the interval [0,T] into 7000 subintervals and integrating using the 
fnint routine. Next, the zeros of Dn{X) ± 2 were calculated by means of the 
fnzeros routine. 

In the following tables, the Mathieu eigenvalues were calculated employing 
the SPPS representation ([T^ for two values of the parameter r. For comparison 
the same eigenvalues from the National Bureau of Standards (NBS) tables are 
also displayed |19l. 





r = 1 


r = 1 


n 


Xn (SPPS) 


A„ (NBS ) 





-0.455139055973837 


-0.45513860 


1 


-0.110248420387377 


-0.11024882 


2 


1.859107160521687 


1.85910807 


3 


3.917024962694820 


3.91702477 


4 


4.371299312651704 


4.37130098 


5 


9.047736927007582 


9.04773926 


6 


9.078369587941564 


9.07836885 


7 


16.033018848985410 


16.03297008 


8 


16.033785039658117 


16.03383234 


9 


25.020598536509114 


25.02084082 


10 


25.021087773318282 


25.02085434 
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r = 5 


r = 5 


n 


A„ (SPPS) 


A„ (NBS ) 





-■5.8()()()4577724278() 


-.").80()()4()()2 


1 


-5.790080596840196 


-5.79008060 


2 


1.858191484309548 


1.85818754 


3 


2.099460384254221 


2.09946045 


4 


7.449142541577460 


7.44910974 


5 


9.236327731534002 




6 


11.548906947651728 




7 


16.648219815375526 




8 


17.096668282587867 




9 


25.510753265631860 


25.51081605 


10 


25.551677357240167 


25.54997175 



Figures 1 and 2 display the plots of the calculated Hill discriminants for two 
values of the Mathieu parameter. 



6 Conclusions 

In summary, in the present work we obtained a representation of the Hill dis- 
criminant in the form of a series in powers of the spectral distance with respect 
to the first eigenvalue and also tested its efficiency for numerical calculations of 
the spectral bands using the Mathieu case. Moreover, we demonstrate that the 
Hill discriminant is invariant under the Darboux transformation generated by 
the first eigenfunction. 
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Figure 1: The polynomial Dn{X) for the Mathieu equation with the parameter 
r = 1 calculated by means of formula i[23|) for N = 100. 
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Figure 2: Same as in the previous figure but for r = 5. The Grst minimum goes 
down to -292.0066. 
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